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@ Frequency-Modulated, Continuous-Wave Laser Ranging Using 
Photon-Counting Detectors 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Optical ranging is a problem of esti- 
mating the round-trip flight time of a 
phase- or amplitude-modulated optical 
beam that reflects off of a target. Fre- 
quency-modulated, continuous-wave 
(FMCW) ranging systems obtain this es- 
timate by performing an interferometric 
measurement between a local fre- 
quency-modulated laser beam and a de- 
layed copy returning from the target. 
The range estimate is formed by mixing 
the target-return field with the local ref- 
erence field on a beamsplitter and de- 
tecting the resultant beat modulation. In 
conventional FMCW ranging, the source 
modulation is linear in instantaneous 
frequency, the reference-arm field has 
many more photons than the target-re- 
turn field, and the time-of-flight esti- 
mate is generated by balanced differ- 
ence-detection of the beamsplitter 
output, followed by a frequency-domain 
peak search. 


This work focused on determining the 
maximum-likelihood (ML) estimation al- 
gorithm when continuous-time photon- 
counting detectors are used. It is 
founded on a rigorous statistical charac- 
terization of the (random) photoelec- 
tron emission times as a function of the 
incident optical field, including the dele- 
terious effects caused by dark current 
and dead time. These statistics enable 
derivation of the Cramer-Rao lower 
bound (CRB) on the accuracy of FMCW 
ranging, and derivation of the ML esti- 
mator, whose performance approaches 
this bound at high photon flux. 

The estimation algorithm was devel- 
oped, and its optimality properties were 
shown in simulation. Experimental data 
show that it performs better than the 
conventional estimation algorithms 
used. The demonstrated improvement is 
a factor of 1.414 over frequency-domain- 
based estimation. 


If the target interrogating photons 
and the local reference field photons 
are costed equally, the optimal alloca- 
tion of photons between these two arms 
is to have them equally distributed. This 
is different than the state of the art, in 
which the local field is stronger than the 
target return. The optimal processing of 
the photocurrent processes at the out- 
puts of the two detectors is to perform 
log-matched filtering followed by a sum- 
mation and peak detection. This implies 
that neither difference detection, nor 
Fourier-domain peak detection, which 
are the staples of the state-of-the-art sys- 
tems, is optimal when a weak local oscil- 
lator is employed. 

This work was done by Baris I. Erkmen of 
Caltech, and Zeb W. Barber and Jason Dahl 
of Montana State University for NASA’s Jet 
Propulsion Laboratory. For more information, 
contact iaoffice@jpl.nasa.gov. NPO-48866 



Several modeling methods are examined. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


For planning of Mars surface mis- 
sions, to be operated on a sol-by-sol 
basis by a team on Earth (where a “sol” 
is a Martian day), activities are de- 
scribed in terms of “sol types” that are 
strung together to build a surface mis- 
sion scenario. Some sol types require 
ground decisions based on a previous 
sol’s results to feed into the activity 
planning (“ground in the loop”), while 
others do not. Due to the differences 
in duration between Earth days and 
Mars sols, for a given Mars local solar 
time, the corresponding Earth time 
“walks” relative to the corresponding 
times on the prior sol/day. In particu- 
lar, even if a communication window 
has a fixed Mars local solar time, the 


Earth time for that window will be ap- 
proximately 40 minutes later each suc- 
ceeding day. Further complexity is 
added for non-Mars synchronous com- 
munication relay assets, and when 
there are multiple control centers in 
different Earth time zones. 

The solution is the development of 
“ops efficiency factors” that reflect the 
efficiency of a given operations config- 
uration (how many and location of 
control centers, types of communica- 
tion windows, synchronous or non-syn- 
chronous nature of relay assets, sol 
types, more-or-less sustainable opera- 
tions schedule choices) against a theo- 
retical “optimal” operations configura- 
tion for the mission being studied. 


These factors are then incorporated 
into scenario models in order to deter- 
mine the surface duration (and there- 
fore minimum spacecraft surface life- 
time) required to fulfill scenario 
objectives. The resulting model is used 
to perform “what-if” analyses for varia- 
tions in scenario objectives. The ops ef- 
ficiency factor is the ratio of the figure 
of merit for a given operations factor to 
the figure of merit for the theoretical 
optimal configuration. 

The current implementation is a pair 
of models in Excel. The first represents 
a ground operations schedule for 500 
sols in each operations configuration for 
the mission being studied (500 sols was 
chosen as being a long enough time to 
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capture variations in relay asset interac- 
tions, Earth/Mars time phasing, and sea- 
sonal variations in holidays) . This model 
is used to estimate the ops efficiency fac- 
tor for each operations configuration. 

The second model in a separate Excel 
spreadsheet is a scenario model, which 
uses the sol types to rack up the total 
number of “scenario sols” for that sce- 
nario (in other words, the ideal number 
of sols it would take to perform the sce- 
nario objectives). Then, the number of 
sols requiring ground in the loop is cal- 


culated based on the soil types con- 
tained in the given scenario. Next, the 
scenario contains a description of what 
sequence of operations configurations is 
used, for how many days each, and this is 
used with the corresponding ops effi- 
ciency factors for each configuration to 
calculate the “ops duration” correspon- 
ding to that scenario. Finally, a margin is 
applied to determine the minimum sur- 
face lifetime required for that scenario. 

Typically, this level of analysis has not 
been performed until much later in the 


mission, and has not been able to influ- 
ence mission design. Further, the notion 
of moving to sustainable operations dur- 
ing Prime Mission — and the effect that 
that move would have on surface mis- 
sion productivity and mission objective 
choices — has not been encountered 
until the most recent rover missions 
(MSL and Mars 2018). 

This work was done by Sharon L. Layback 
of Caltech for NASA’s fet Propulsion Labora- 
tory. Further' information is contained in a 
TSP (see page 1 ). NPO-48262 


Q GPU Lossless Hyperspectral Data Compression System 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Hyperspectral imaging systems on- 
board aircraft or spacecraft can acquire 
large amounts of data, putting a strain 
on limited downlink and storage re- 
sources. Onboard data compression can 
mitigate this problem but may require a 
system capable of a high throughput. In 
order to achieve a high throughput with 
a software compressor, a graphics pro- 
cessing unit (GPU) implementation of a 
compressor was developed targeting the 
current state-of-the-art GPUs from 
NVIDIA®. 

The implementation is based on the 
fast lossless (FL) compression algorithm 
reported in “Fast Lossless Compression 


of Multispectral-Image Data” (NPO- 
42517), NASA Tech Briefs, Vol. 30, No. 8 
(August 2006) , page 26, which operates 
on hyperspectral data and achieves ex- 
cellent compression performance while 
having low complexity. The FL compres- 
sor uses an adaptive filtering method 
and achieves state-of-the-art perform- 
ance in both compression effectiveness 
and low complexity. The new Consulta- 
tive Committee for Space Data Systems 
(CCSDS) Standard for Lossless Multi- 
spectral & Hyperspectral image com- 
pression (CCSDS 123) is based on the 
FL compressor. The software makes use 
of the highly-parallel processing capa- 


bility of GPUs to achieve a throughput 
at least six times higher than that of a 
software implementation running on a 
single-core CPU. This implementation 
provides a practical real-time solution 
for compression of data from airborne 
hyperspectral instruments. 

This work was done by Nazeeh I. Aranki, 
Didier Keymeulen, Aaron B. Kiely, and 
Matthew A. Klimesh of Caltech for NASA’s Jet 
Propulsion Laboratory. For more information, 
contact iaoffice@jpl.nasa.gov. 

The software used in this innovation is avail- 
able for commercial licensing. Please contact Da n 
Broderick at Daniel.F.Broderick@jpl, nasa.gov. 
Refer to NPO-48571. 


@ Robust, Optimal Subsonic Airfoil Shapes 

Ames Research Center, Moffett Field, California 


A method has been developed to cre- 
ate an airfoil robust enough to operate 
satisfactorily in different environments. 
This method determines a robust, opti- 
mal, subsonic airfoil shape, beginning 
with an arbitrary initial airfoil shape, 
and imposes the necessary constraints 
on the design. Also, this method is flex- 
ible and extendible to a larger class of 
requirements and changes in con- 
straints imposed. 

In one embodiment, process steps in- 
clude providing a specification of a de- 
sired pressure value at each of a se- 
quence of selected locations on the 
surface of a turbine airfoil; providing an 
initial airfoil shape; providing a state- 
ment of at least one constraint to which a 
final airfoil shape must conform; using 
computational fluid dynamics (CFD) to 
estimate a pressure value at each of the 


selected perimeter locations for the ini- 
tial airfoil shape; using CFD to deter- 
mine the pressure distribution for the 
airfoil shapes that are small perturba- 
tions to the initial airfoil shape; and 
using an estimation method, such as a 
neural network, a support vector ma- 
chine, or a combination thereof, to con- 
struct a response surface that models the 
pressure distribution as a function of the 
airfoil shape using the CFD data. Other 
process steps include using an optimiza- 
tion algorithm to search the response 
surface for the airfoil shape having the 
required pressure distribution, and pro- 
viding at least one of an alphanumeric 
description and a graphical description 
of the modified airfoil shape. 

Constraints may be drawn from the 
following group, or may be one or more 
other suitable constraints: vortex shed- 


ding strength from the trailing edge of 
the airfoil is no greater than a selected 
threshold value; a difference between 
any resonant frequency of the airfoil 
and the vortex shedding frequency is at 
least equal to a threshold frequency dif- 
ference; mass of the airfoil is no larger 
than a threshold mass value; and pres- 
sure value at each of a sequence of se- 
lected locations along the surface of the 
airfoil differs from a corresponding ref- 
erence pressure value by no more than a 
threshold pressure difference value. 

This work was done by Man Mohan Rai of 
Ames Research Center. Further information is 
contained in a TSP (see page 1 ). 

Inquiries concerning rights for the commer- 
cial. use of this invention should be addressed 
to the Ames Technology Partnerships Division 
at 1-855-NASA-BIZ (1-855-6272-249). Refer 
to ARC- 145 86-2 
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